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We investigate theoretically the effects of surface states over the magnetoresistance of Ni-based 
organic spin-valves. In particular we perform ab initio electronic transport calculations for a benzene- 
thiolate molecule chemically attached to a Ni [001] surface and contacted either by Te to another Ni 
[001] surface, or terminated by a thiol group and probed by a Ni STM tip. In the case of S- and Te- 
bonded molecules we find a large asymmetry in the spin-currents as a function of the bias, although 
the I-V is rather symmetric. This leads to a smooth although not monotonic dependence of the 
magnetoresistance over the bias. In contrast, in the case of a STM-type geometry we demonstrate 
that the spin-current and the magnetoresistance can be drastically changed with bias. This is the 
result of a resonance between a spin-polarized surface state of the substrate and the d-shell band-edge 
of the tip. 

PACS numbers: 



I. INTRODUCTION 

The electron spin has recently made its appearance in 
organic electronics, generating a new field which takes 
the suggestive name of molecular-spintronics The 
main idea is to investigate how spin-currents propagate 
through organic molecules possibly affecting their inter- 
nal state. Example of such spin-molecular devices are 
spin-valves made from carbon nanotubes 2], organic in- 
sulators 3] or conducting polymers [4|, |5j , fempto-second 
spectroscopy of molecular bridges connecting semicon- 
ductor nanodots |(| and three terminal devices using 
magnetic molecules as active elements 0, • 

Organic molecules present genuine advantages over 
conventional inorganic metals and semiconductors. Typ- 
ically the spin-orbit interaction, which scales with the 
atomic number Z as Z 4 , is weak in molecules. For in- 
stance the spin-orbit splitting of the valence band of car- 
bon in its diamond form is only 13 meV [jj, to be com- 
pared with 340 meV of GaAs. This yields rather long 
spin-diffusion lengths (l s ) with conservative estimates of 
130 nm in multiwall carbon nanotubes 0, 200 nm in 
polymers and 5 nm in Alq3 monolayers |10|. It 
is important to note that these values are usually ex- 
tracted from the magnetoresistance (MR) of spin-valves 
using the Jullier's formula ^lj an( i therefore also in- 
clude possible spin-flip at the surface of the magnetic 
electrodes. This means that the intrinsic spin-diffusion 
length is such organics may be considerably longer. Fur- 
thermore one should remember that at least in polymers 
and Alq3 the resistivity is rather large and therefore the 
spin-relaxation times associated to the values of l s given 
above are quite long. 

In organic molecules also the hyperfine interaction, an- 
other possible source of spin-decoherence, is weak. Most 
of the molecules used for spin-transport are 7r-conjugate 
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and the conduction is mostly through molecular states 
with amplitude over the carbon atoms. Carbon, in its 
most abundant isotopic form, 12 C, has no nuclear spin, 
and therefore is not hyperfine active. Moreover the ir- 
states are usually delocalized and the interaction is in 
general weak even over those nuclei with non-vanishing 
nuclear spin. 

In addition to the long spin-diffusion length organic 
materials offer a range of new physical properties, diffi- 
cult to find and tune in conventional inorganic materials. 
For instance the electrical conductivity of a polymer can 
be changed by over twelve orders of magnitude with the 
incorporation of various counterions |12j . More inter- 
estingly the current-carrying elementary excitations are 
not band-like but usually are in the form of electronic- 
vibronic correlated states, with peculiar spin character- 
istics [T3| . Furthermore magnetic molecules with collec- 
tive high-spin ground states are available 0|. In these 
the transport can be dominated by long-living metastable 
states, giving rise to spin-dependent non-linearities in the 
I-V characteristics 00. 

Finally a less explored avenue is that of exploiting 
the many alternatives for bonding a molecule to a mag- 
netic surface. The electronic structure of the anchoring 
group can in fact change quite drastically the magneto- 
transport response of a molecule [l5| paving the way for a 
chemical design of magnetic devices. This is precisely the 
problem we address in the present work. In particular we 
investigate how spin-polarized surface states, almost al- 
ways present in the case of organic spin- valves, can affect 
the transport. 



II. COMPUTATIONAL METHOD 

All the calculations presented in this work have been 
performed with the code Smeagol [TH ll»| . Smeagol 
combines the non-equilibrium Green's function (NEGF) 
method 0,^3 with the Kohn-Sham 0] form of density 
functional theory (DFT) [20|. In its present implemen- 
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tation Smeagol uses the code SIESTA [21| as the main 
DFT platform. This has the advantage of expanding 
the Kohn-Sham problem over an efficient linear combi- 
nation of localized numerical orbitals [2^, [2j|, making 
large systems tractable without the need of massive com- 
putational resources. In addition core electrons are re- 
moved from the self-consistent calculation and described 
by standard norm-conserving scalar relativistic Troullier- 
Martins pseudopotentials. 

The central quantity that Smeagol evaluates is the non- 
equilibrium Green's function G(E) for the portion of the 
device where the potential drops. Usually this is called 
"extended molecule" (EM) and comprises the molecule 
itself and part of the current /voltage probes. Such a 
Green's function simply reads 
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In our simulations the EM includes three Ni planes at 
each side of the molecule and the cross section is large 
enough to prevent intermolecular interaction. We use pe- 
riodic boundary conditions in the direction orthogonal to 
the transport with a 9 k-points in the two-dimensional 
Brillouin zone. The real-space grid cut-off is fixed at 
450 Ry [2l| and multiple-C basis sets are constructed in 
the following way: double-C for the s shell of H, C, S, 
Te and Ni, and for the Ni d shell, double-C plus polar- 
ization for the p shells of C, S and Te, and single-^ for 
Ni p orbitals. Finally the charge density is obtained by 
integrating the Green's function over 96 imaginary and 
300 real energy points . 



III. RESONANT TRANSPORT IN ORGANIC 
SPIN-VALVES 



where Hpm is the Hamiltonian of the EM, E is the en- 
ergy and El (Er) is the self-energy for the left-hand side 
(right-hand side) electrode. These latter ones are non- 
hermitian energy-dependent matrices describing the in- 
teraction of the EM with the current/ voltage electrodes 
and effectively allow us to map a hermitian problem 
for an open non-periodic system onto a non-hermitian 
problem for a finite system. They can be written as 
s « = Hl EM g a H aEM with H aEM the coupling matrix 
between the leads a (a —L, R) and the EM. The g a are 
surface Green's functions for the leads, i.e. the Green's 
function for a semi-infinite lead evaluated at the termina- 
tion plane. They are constructed using a semi-analytical 
expression [25|. after removing the singularities of the 
coupling matrices with a renormalization scheme fl5| . 

Importantly when Hpm is a function of the charge den- 
sity p, as in DFT, then a simple prescription for a self- 
consistent calculation can be designed. In fact, a close 
formula for the non-equilibrium charge density in terms 
of the Green's function G(E) is available [T{|, and can 
be used in conjunction with equation for deriving the 
electronic structure of the EM in presence of the leads. 
External bias can be added in the self-consistent proce- 
dure by imposing the appropriate boundary conditions to 
the Hartree potential in the leads. This simply produces 
a rigid shift of the whole spectra in the current /voltage 
electrodes, since they are good metals and therefore local 
charge neutral. 

Finally the two probe current can be extracted by using 
a Landauer-likc formula 
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where G = G(E), T a = i[E a - EjJ, f a (x) is the Fermi 
function and /xl/r = Ep ± eV/2 with Ep the Fermi en- 
ergy of the leads and V the external bias. Note that ulti- 
mately the current is nothing but the integral of the bias- 
dependent transmission coefficient T(E) — GTlG^Tr 
evaluated over the bias window. 



Consider a typical spin- valve constructed from a strong 
ferromagnet and a non-magnetic molecule. In this set- 
up the transport is determined by the alignment of the 
molecular levels with the Fermi energy (Ep) of the elec- 
trodes 0, 0| . Let us consider the situation in which 
the transmission through the molecular level is small, 
but the anchoring groups form a spin-polarized localized 
state with the d-orbitals of the magnetic electrodes. The 
spin-splitting of such state will be comparable with the 
band-exchange of the ferromagnet, although its actual 
form and position will depend on the details of the an- 
choring structure. 
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FIG. 1: Cartoon representing the alignment of the spin- 
polarized surface states associated to the two magnetic elec- 
trodes forming a spin-valve: (a) parallel alignment and for- 
ward bias, (b) antiparallel alignment and forward bias, (c) 
antiparallel alignment and reverse bias. Panels to the left are 
for V—0 and the bias increases from left to right. The red 
boxes mark the three resonant conditions. 

In this condition of strong coupling the position of such 
surface state (SS) is strongly pinned at the chemical po- 
tential of the relevant electrode and shifts rigidly with 
the bias. Then, the alignment of the SS associated to 
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the left electrode of the device with that of the right 
determines the current. In particular for identical elec- 
trodes and anchoring groups there are three different res- 
onant conditions (see figure^) depending on the bias and 
the mutual alignment of the magnetization vectors in the 
contacts. In the case of parallel alignment (PA) the state 
associated to the left-hand side electrode is lined up with 
that of the right only at zero bias and for both spin di- 
rections (left panel of figure^,). If it is also located at 
Ep, then the current will be large around V=0. A finite 
bias reduces the overlap between the states, regardless of 
the bias polarity (in the picture we report forward bias 
only). This means that the polarization of the currents 
decreases with bias and approximately vanishes for a bias 
voltage as large as the state width. 

The other two resonant conditions are found for the 
antiparallel alignment (AA) of the electrodes. In fact 
the SSs on the two sides of the device overlap for a bias 
of the order of the exchange splitting of the state itself. 
For forward bias such overlap is between the majority 
component of the SS to the left with the minority of that 
to the right (last panel to the right of figure \T]p) . The 
opposite is for reverse bias (last panel to the right of 
figure [1J;). This gives the interesting result that for the 
AA both the spin-currents are highly asymmetric as a 
function of bias, although the total current is not. Thus, 
if we define the MR ratio as Rmr = {Ipa — Xaa)/-Taa, we 
will expect a minimum in Rmr for biases of the order of 
the SS exchange splitting. This occurs in correspondence 
of a large current in the AA configuration. 

A completely different situation is encountered when 
the attachment of the molecule to a surface is highly 
asymmetric. In this case different resonance conditions 
will be present depending on the details of the struc- 
ture and both the spin-current and the MR might have 
non-trivial bias dependence |26|. Similar arguments 
have been also used for explaining rectification in non- 
magnetic molecular junctions using electrodes with d 
states close to Ep [27| . 



IV. S-BENZENE-Te BETWEEN Ni SURFACES 

We start our analysis by investigating a benzene 
molecule chemically attached to the hollow site of two Ni 
[001] surfaces. As anchoring groups we use S on one side 
and Te on the other. These both form a stable chem- 
ical bond with Ni, although the two bond lengths are 
quite different (1.28 A for S and 1.77 A for Te) reflecting 
the rather different atomic radii of the anchoring atoms. 
With this configuration we expect an asymmetric I-V 
characteristics, although we do not expect any rectifica- 
tion, since the coupling with the electrodes is strong. 

In figure |21 we present the I-V characteristics for both 
the PA and AA, the MR ratio as a function of the bias 
voltage and the transmission coefficients as a function of 
energy for zero bias in the PA configuration. From the 
transmission coefficient it is clear that most of the cur- 
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FIG. 2: Transport properties of a S-Benzene-Te molecule 
sandwiched between Ni [001] surfaces. In the large panels the 
I-V characteristic for both the PA (a) and A A (b). In the left 
inset the MR ratio as a function of bias, and in the right one 
the transmission coefficient as a function of energy for the PA 
and zero bias. 



rent through the device is carried by one broad transmis- 
sion peak located just below the Fermi level (Ep) for the 
majority spins and just above for the minority. These 
are associated with a hybrid SS formed from Ni-d or- 
bitals with the p-shells of both S and Te. The resulting 
transport properties are those expected from the simple 
discussion of the previous section. 

In the PA the two spin-currents are rather similar 
since no resonant condition is ever met except for zero 
bias. Note also that both spin-currents almost satu- 
rate at around V = ±0.5 Volt, which approximately 
corresponds to the width of the SS. Finally it is clear 
that the spin polarization of the current, defined as 
P = (JT - I l )/(P + I 1 ) (I a is the spin a current), de- 
creases as the bias increases and approximately vanishes 
for V -1 Volt. 

In the AA the two spin-currents individually present 
a rather asymmetric bias dependence reflecting the two 
resonant conditions. However the total current is almost 
completely symmetric and the tiny asymmetry is asso- 
ciated with the slightly different bonding of S and Te 
with the Ni hollow site. Importantly also in this case the 
current saturates at around V = ±0.5 Volt, when the 
resonant condition is not met any longer. This gives a 
MR ratio with a minimum and negative sign at around 
V = ±0.5 Volt. Interestingly the spin-resolved MR ratio 

r mr = ( J pa- J Aa)/ J AA (note that R MR + R MR = Rmr) 
is rather asymmetric with V, reflecting the asymmetry of 
the spin-currents. Note finally that the MR ratio changes 
sign when the bias is increased from zero, since we gradu- 
ally move from the resonant condition for the PA to those 
for the AA. Similar results have been recently reported 
for different anchoring structures [2^ and molecules [2^] . 
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V. STM CONFIGURATION 

The junction described in the previous section presents 
rather large currents and the typical transmission coeffi- 
cient approaches T = 1. This is not the ideal situation 
for studying resonances, which in contrast dominate the 
transport in the case of overall small transmission. For 
this reason we move to investigating a STM-type geom- 
etry. Here we consider the same benzene molecule at- 
tached to the hollow site of a Ni [001] surface via S and 
terminated with a thiol group. This latter is then in 
loose contact with a single-atom Ni STM tip positioned 
on top of S (with a Ni-S separation of 4 A) and shifted 
perpendicularly to the benzene plane by ml. 2k. (see inset 
of figure In this configuration we probe the p z orbital 
of S (i.e. the one perpendicular to the benzene plane), 
which is involved in the transport. 
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FIG. 3: I-V characteristic for a benzene-thiol-thiolate 
molecule attached to the hollow site of Ni [001] and probed 
with a Ni STM tip. In the insets the geometrical configuration 
of the device and Rmr as a function of bias. 

In figure|31we present both the I-V characteristics and 
the MR for such a device. At variance with the previous 
case the total currents for the PA and AA are rather 
similar for forward bias, but differ substantially for the 
negative voltages. In this case the PA current drastically 
increases for V between -0.1 Volt and -0.5 Volt and then 
saturates. In turns this produces a large positive peak 
in the MR reaching 70% for V around -0.5 Volt. Away 
from this peak the MR ratio is a smooth function of the 
bias and generally it is quite small. 

In order to have a better understanding of this second 
situation in figure 0] we show the orbital resolved density 
of states for our device, where in particular we project 
over the Ni [001] surface, the benzene molecule and the 
d orbitals of the tip. From the picture we can identify 
the two spin-polarized states mainly responsible for the 
transport. The first is the tip d-orbitals band-edge (red 
circles in figure 0J) with majority and minority compo- 
nents approximately 200 meV below the corresponding d 
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FIG. 4: Orbital resolved density of state for benzene- 
thilolate-thiol molecule attached to Ni [001] and probed by 
a Ni STM tip. The tip and surface states are indicated re- 
spectively with red circles and blue squares. 

with spin-polarized components (blue squares) approxi- 
mately 200 meV above the corresponding d band-edge of 
the substrate. 

It is then clear that a resonant condition is found when- 
ever the tip d band-edge is aligned with SS of the sub- 
strate. However, in contrast to the cartoon of figure ^ 
the resonant states of the left-hand side electrode (asso- 
ciated to the Ni [001] substrate) and of right-hand side 
electrode (associated to the tip) are in different positions. 
The alignment that we obtain from figure ^ immediately 
suggests that, regardless of the magnetization alignment 
of the electrodes, there cannot be any resonant condition 
for positive bias (current going from the substrate to the 
tip), since the tip state is lower than the SS of the sub- 
strate. In contrast a resonant condition can be matched 
at negative bias for PA and both spins. 

In figure[3]we demonstrate this resonant condition and 
present the transmission coefficient as a function of en- 
ergy in the PA for various biases. We note that for zero 
bias T(E) displays two 0.5 eV wide peaks around Ep 
for both spin directions. These directly reflect the tip 
d orbital band-edge and the spin-polarized SS. As the 
negative bias increases the tip and substrate state move 
to resonance as demonstrated by the large transmission 
peak at V ~ -0.4 Volt. This occurs for both majority 
and minority spins, however only the minority contribute 
to the current since only the minority peak falls within 
the bias window. Finally as the bias keeps increasing 
the resonance is destroyed and the enlargement of the 
bias window is compensated by a drastic reduction of the 
transmission peaks. This causes the current to saturate. 

Note that for positive bias no resonant behavior is ob- 
served within the bias window at any of the voltages in- 
vestigated, although some resonances are observed ap- 
proximately 1 eV below E-p. A similar non-resonant sit- 
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FIG. 5: Transmission coefficient as a function of energy at 
different bias for benzene-thilolate-thiol molecule attached to 
Ni [001] and probed by a Ni STM tip: parallel alignment only. 
Note the sharp increase of the transmission for V~ -0.4 Volt. 
The solid (dashed) line represent respectively majority and 
minority spins. The vertical red lines indicate the bias win- 
dow. 



uation is encountered for the AA case and it is not pre- 
sented here. 



VI. CONCLUSIONS 

In conclusion we have investigated resonant transport 
in organic spin-valves. In the case of a single resonant 
state we predict a maximum in the current for the an- 
tiparallel alignment of the electrodes at a bias compa- 
rable with the exchange-splitting of the resonant state. 
In contrast in the case of highly asymmetric anchor- 
ing structures multiple resonant states may appear and 
different resonant conditions can be found. In particu- 
lar we investigate the almost symmetric case of benzene 
molecule chemically bonded to Ni [001] surface via S and 
Te, and that of a STM-type geometry. Importantly the 
MR ratio is a non-trivial function of bias and can as- 
sume both positive and negative values. Therefore such 
a strong variation of Rmk with bias seems to be a dis- 
tinctive feature of organic spin- valves. 
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